Characterization of 2D fermionic insulating states 
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Inspired by the duality picture between superconductivity (SC) and insulator in two spatial di- 
mension (2D), we conjecture that the order parameter, suitable for characterizing 2D fermionic 
insulating state, is the disorder operator, usually known in the context of statistical transformation. 
Namely, the change of the phase of the disorder operator along a closed loop measures the particle 
density accommodating inside this loop. Thus, identifying this (doped) particle density with the dual 
counterpart of the magnetic induction in 2D SC, we can naturally introduce the disorder operator 
as the dual order parameter of 2D insulators. The disorder operator has a branch cut emitting from 
this "vortex" to the single infinitely far point. To test this conjecture against an arbitrary 2D lattice 
models, we have chosen this branch cut to be compatible with the periodic boundary condition 
and obtain a general form of its expectation value for non-interacting metal/insulator wavefunction, 
including gapped mean-field order wavefunction. Based on this expression, we observed analytically 
that it indeed vanishes for a wide class of band metals in the thermodynamic limit. On the other 
hand, it takes a finite value in insulating states, which is quantified by the localization length or 
the real-valued gauge invariant 2-from dubbed as the quantum metric tensor. When successively 
applied along a closed loop, our disorder operator plays role of twisting the boundary condition of 
a periodic system. We argue this point, by highlighting the Aharonov-Bohm phase associated with 
this non-local operator. 

PACS numbers: 74.81.-g, 71.10.Pm, 71. 23. An, 77.22. Ch, 72.80.Sk 



I. INTRODUCTION 



The high-T c cuprates and the colossal magnetoresis- 
tive manganites, the two main realms in condensed mat- 
ter physics, are both doped correlated insulators. As 
type II superconductors accommodate magnetic flux by 
allowing a spatially inhomogeneous distribution of the 
superconducting order parameter, such fermionic insu- 
lating states do not necessarily remain spatially uniform 
against doping. Indeed recent STM experiments have 
observed that a hole rich region in a doped Mott insu- 
lator such as Bi-2212 forms a nanoscale granular island 
on top of the insulating background^. A similar kind 
of nanoscale electronic inhomogeneity was also found in 
Ca2- :c Na a ;Cu02Cl2, where doped holes form a checker- 
board pattern 2 .. On the other hand, extensive Lorentz 
optical microscopy measurements on doped (e.g., mixed 
valence) manganites have revealed that charge ordered 
and ferromagnetic metallic patches are phase-separated 
in such relatively large length scales as a few microme- 
ter—. The question naturally arises as follows : Is there 
some universal classification of insulators in terms of their 
different types of behavior against doping? The other 
way around, if it exists, what kind of microscopic ingre- 
dients would determine these differences? 

These experimental questions/observations as well as 
almost all the other experimental observables are di- 
rectly accessible through the "local electronic polariza- 
tion" P(r), which acquires finite mass only in dielectrics. 
Contrary to other observables, however, the local elec- 
tronic polarizations cannot be uniquely determined by 



those vector fields which are divergence-free. Namely, 
using arbitrary analytic scalar functions 0(f), we can in- 
troduce those vector fields which have nothing to do with 
the local charge density p(r) = V • P(r); 



P'(f) = P(r) + xV0(f), 
V • P'(r) = V • P(r). 



(1) 



Due to this arbitrariness, constructing effective theories 
of dielectrics in such a way that the local polarization 
becomes explicit has been one of the long-standing issues 
in condensed matter theories community, while widely 
demanded from experimental sides on a general ground. 

The central idea of the duality picture is to regard this 
local electronic polarization as the "gauge field" which 
is intrinsic in matters and to treat its divergence-free 
part as the unphysical "gauge" degree of freedomi^i^ 
To be more specific, we suppose that there exist order 
parameters, implicit in a microscopic model, which are 
coupled with the local electronic polarization such that 
its condensation makes this "gauge fields" massive via 
the Higgs mechanismi^. Accordingly, the "gauge" in 
eq.([T]) and that of this order parameter 77(f) are specified 
in a set, as in superconductors; 



i]'(r) = 77(f) e 



i9(r) 



(2) 



where the vortex-free scalar function 0(r) should be iden- 
tical to that in Eq. (JTJ. 

In this paper, observing this duality picture, we will 
introduce this complex-valued order parameter 77(f) ex- 
plicitly in terms of original fermion ( electron ) operators. 
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Then we calculate the expectation value of this order pa- 
rameter with respect to some simple wave functions, so 
that we can argue that this order parameter has indeed 
a finite amplitude, i.e. condensed, only in fermionic in- 
sulators. 

This paper is organized as follows : In Sec. II, we look 
further into the duality between superconductivity and 
insulator only to arrive at the conjecture that the appro- 
priate order parameter for a 2D fermionic insulating state 
is the disorder operator (DOP) defined as eqs. (|4l9p . We 
then give a general expression to the expectation value 
of the DOP. In Sec. Ill, we will see that this expecta- 
tion value indeed vanishes for band metals having vari- 
ous kinds of Fermi surface (F.S.) in the thermodynamic 
limit. In Sec. IV, we will consider the opposite limit, 
i.e., the case of a band insulator/gapped mean-field or- 
der state close to the atomic limit. Thereby, we observe 
that the expectation value of the DOP is in turn char- 
acterized by the localization length and thus remain fi- 
nite. Then these two observations, i.e. those in sec. Ill 
and in sec. IV, lead us to extrapolate the behaviour of 
the DOP in general insulating states. Sec. V is devoted 
to the arguments on the relation between the DOP and 
"twisting boundary condition" , the latter known to give 
a definite criterion for insulating states in arbitrary spa- 
tial dimensio n 10 ' 11 ! 12 ' 20 . Sec. VI contains not only the 
brief summary (see Table. II) but also discussions on the 
behaviour of our DOP in the off-diagonal long ranged or- 
dered states. We also mentioned there about the possible 
microscopic candidates of the counterpart of the magnetic 
penetration depth/ coherence length (see Table. I) and the 
future possible progress. Some details of the calculation 
are left to the appendix. 



II. THE DISORDER OPERATOR AND THE 
INSULATING ORDER PARAMETER 



On duality between superconductivity and 
insulator 



area S. Then, we have, 



$ <t> V(f>(r) ■ dl = I b(r) d 2 r, 



i)S 



Pes 



(3) 



where V = (d/dx, d/dy) and dS is the boundary of this 
area S. We now invoke this relation to identify the dual 
quantity of A(f). 



2D electronic insulator 


2D superconductor 


doping 
doped particle density: 
Sp(r) = p(r) - p 


applied magnetic field 
magnetic induction: 
6(f) 


"type I insulator (?)" 
"type II insulator (?)" 

? 

? 


type I superconductor 
type II superconductor 
penetration depth 
coherence length 


disorder operator: 
n(r) = |j7(r)|e ie(F) 


order parameter: 
A(f) = |A(f)|e^ (?) 


a„(f) =e^P v (r) = V M 0(f) 
8p(r) = V • P = V x a 


A,(r) = V M 0(f) 
56(f) = V x A 



TABLE I: The duality between 2D insulator and supercon- 
ductivity 

In the duality relation between 2D superconductivity 
and 2D insulator, applying magnetic field in supercon- 
ductors corresponds to doping in insulators 4 (Table. I). 
Thus the doped particle density 5p(f) = p(r) — p corre- 
sponds to magnetic flux in superconductors. Physically 
speaking, 5p(r) is obtained from a local electronic polar- 
ization P(r); 5p(r) = V ■ P(r). Then introducing the 
"gauge field" a(f) such that a M (f) = e IJ , 1/ P v (r), we have 
only to find out a scalar function whose gradient is this 
"gauge field", i.e. d / j,6(r) = a M (f). Namely, we can iden- 
tify this scalar function 9(f) as the phase part of the dual 
order parameter, i.e. counterpart of </>(f). We thus reach 
the conjecture that the disorder operator (DOP) fj(r)^ 
should play role of an insulating order parameter; 



Let us denote the usual superconducting order param- 
eter as A(f) = |A(r)|e l< ^ r ). In this paper, we consider 
only two spatial dimension (2D), i.e. r — (x,y). If the 
magnetic penetration depth A is large enough compared 
with the coherence length £ (if n = A/£ > l/y/2 in the 
conventional Ginzburg-Landau theory^), then the sys- 
tem shows a type II superconducting behavior, accom- 
modating magnetic flux b(r) pinned to vortices in the 
system. The amplitude of order parameter |A(f)| is spa- 
tially not uniform and can have zeros. Then its phase 
<j)(f) acquires an ambiguity of integer multiples of 2-zr 
around its zeros. The holonomy of this phase around 
such a vortex is proportional to the number of magnetic 
flux pinned to this vortex, counted in units of the flux 
quantum $o = hc/(2e). Let a flux penetrate a specific 



fj(z) = exp 



dVlog(z' -z){p(z l )-p} 



z'^z 



(4) 



Here we have introduced a complex variable z = x + iy, 
and d 2 z should be understood to be d 2 z — dzdz/(2i). 
Note that fj(r) is defined on a 2D infinite plane without 
any boundary condition. Using this explicit form, one 
can readily verify that eq. |(4]) indeed satisfies our leading 
principle, 

-!- / Varg 77(f) • dl= { {p(r) - p} d 2 f. (5) 
2lT JdS J Pes 

This states that the winding number associated with the 
dual order parameter 77(f) is identical to the doped par- 
ticle number inhabiting within 5*. 



3 



B. Disorder operator for a periodic lattice 

Table. I shows the correspondence between 2D insu- 
lator and superconductivity. In the middle two empty 
seats are reserved for the unknown counterparts of the 
magnetic penetration depth and the coherence length in 
2D insulator. If we push forward with this duality rela- 
tion, our eventual goal might be to construct a Ginzburg- 
Landau theory for 2D (doped) insulating state in terms 
of this non-local operator, which unambiguously lets us 
complete this table quantitatively. As a first step to- 
ward this direction, we will demonstrate in the following 
sections that the expectation value of the DOP 

1. indeed vanishes in various types of 2D band met- 
als, i.e., the DOP is not condensed (see Sec. Ill), 
whereas 

2. it remains finite in 2D band insulators, i.e., the 
DOP is condensed (see Sec. IV). 

If one tries to make sure these statements, one might 
immediately notice that the DOP defined in eq. ((4]) in 
itself is incompatible with the periodic boundary condi- 
tion (PBC) which we usually presume. The branch cut 
of logarithm in eq. has only a single end point, while, 
mathematically speaking, it is impossible to embed such 
a branch cut into a torus. 

A possible and the most plausible way out would be 
to replace z' — z in this logarithm by a doubly periodic 
function CO' —x,y' — y), i.e. CO) V + L) = CO + L,y) = 
CO, y), which reproduce z' — z when \z' — z\ <C L, 

lim CO' - x,y' - y) ~ z' - z. (6) 

\z' — z\/L— >0 

Here L denotes the linear dimension of a system size. 
As the simplest function satisfying this requirement, we 
consider in this paper the following function, 

CO' - x,y' -y) = 
-i(e i2 r (x ' ~ x ~ Xo) - l) + (e^v'-v-wl - 1). (7) 

Apart from the continuum variables z and z' used in 
eq. ([H, r = (x, y) and r* = 0', y') in eq. (JT]) are defined to 
take a discrete value, specifying lattice points in a partic- 
ular lattice model, (xo, yo) was then introduced, so as to 
avoid the logarithmic singularity when log C are summed 
with respect to r 1 over all lattice points (see eq. (0) ). In 
a simple square lattice with its lattice constant "a" , we 
will take (xo,yo) = (§, §) (see Fig. [T|). This function, i.e. 
eq. ([7]), in fact satisfies the requirement eq. ([6]). Namely, 
when \z' — z\ <C L, it actually reads, 

CO' - x, y' - y) ^ -j-{0' -x-x ) + i(y' -y- y )}- 

It is also a doubly periodic function. As a result, when 
seen as a continuous function of (x 1 , y'), CO' ~ x > v' — V) 
has two zeros: 

(x',y') = (x + x ,y + y ), (x + x + j,y + y - ^)-( 8 ) 
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FIG. 1: (color online) 2D lattice models. (x,y) and (x',y') 
specify a lattice point (filled circle mark). A unit cell is de- 
picted as bold line in every figure and we use different colors 
(shapes) in order to distinguish Nb number of inequivalent 
lattice points within it. Crossed mark represents a dual lat- 
tice site on which (x + xo,y + yo) is located. (From top to 
down): Triangle and Square lattice with Nb — 1, honeycomb 
and checkerboard lattice with Nb = 2 and kagome lattice 
with Nb = 3. Note that (x + xq, y + yo) (crossed mark) never 
coincides with the original lattice sites on which a fermion 
operator is defined. 

The first zero corresponds to the original vortex intro- 
duced in eq. (g]). The latter one turns out to be the 
antivortex, which was supposed to be located at the sin- 
gle infinitely far point in eq. ([4J. To be more specific, 
when (a/, y') moves around the former (latter) zero point 
anti-clockwise, then logC picks up a phase 2-7T (— 2tt), as 
shown in Fig.O Accordingly, logC has a branch cut run- 
ning from the vortex to the antivortex as in Fig. [3] In 
general, the singularities must appear in pairs in a sys- 
tem obeying PBC. In order to study the DOP explicitly 
in non- interacting states and mean-field ordered states, 
it is necessary to write down a ground state wavefunc- 
tion explicitly. Without PBC, we would not be able to 
do this. We thus conjecture in the remainder of the pa- 
per that the following operator with C defined in eq. ([7]) 
is the appropriate form of the disorder operator (DOP) 
compatible with PBC, 

rj(r) = exp [ ^ log CO' - x,y' - y){p(P) - p}] . (9) 

x',y' 
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(7t/2,-Jt/2) 



FIG. 3: : A branch cut runs from the vortex ((x,y) — (0,0)) 
to the anti-vortex ((x,y) = (f,-f)). Arglog[-i(e" - 1) + 
(e ly — 1)] takes — ^— in the shaded side of this branch cut, 
while takes — ^+ in the other side. 



The general expression for expectation values of 
the disorder parameter 



FIG. 2: (color online) Arglog[-j(e" - 1) + (e iy - 1)] as a 
function of x and y, where we only show the principle branch 
which ranges from — ^ to — ^ in this picture and those below. 



Notice that the ^'-integral in eq. (|4|) avoids z' = z, 
while the summation w.r.t. r 1 = (x',j/) in eq. ([9|) is 
taken over entire lattice points without any restriction. 
As we mentioned above in the square lattice case with its 
lattice constant "a" , this simplification becomes possible, 
just because we have introduced (x ,yo) — (§, §) so that 
log(,(n x &, n y a) is always finite for arbitrary integer n x 
and n y . All 2D lattice models, however, can be also re- 
garded as having a square unit cell (See Fig. [IJ. Then, 
we can always choose (xq, yo) appropriately, such that the 
vortex and antivortex introduced in eq. f5P never coincide 
with any lattice points on which p(r*) is defined (see some 
examples shown in Fig. Qp. Provided that log^r* — f) in 
eq. §§§ takes finite value for an arbitrary lattice point r*, 
then the following arguments do not depend seriously on 
the specific choice of (xo, yo)- 

Before giving the general expression to the expectation 
value of our DOP defined in eqs. (|7l9p . let us fix our 
ideogram of this paper now. Firstly, we always refer a 
lattice constant of a unit cell (square box depicted by 
bold line in Fig. [T|) to "a" . We take a number of this 
unit cell in an entire system to be N x N and thus the 
linear dimension of a system size is given by L — Na. We 
call as Nb, the total number of Bloch bands/inequivalent 
sites within each unit cell, while total number of electrons 
which a ground state w.f. has is denoted by N e . Then the 
average particle number per site p introduced in eq. © 
is given as follows, 



Nb ■ N 2 



(10) 



For a non-interacting system and arbitrary mean-field 
ordered state with N e fermions, a ground wavefunction 
is obtained just by filling the lowest N e one-body states 
which we name simply as a = 1, 2, • • • , N e , 



|GS) 



c[\Q). 



(11) 



Here |0) denotes the fermion-vacuum state. To evaluate 
our DOP for such a ground state, let us apply our DOP 
onto a creation or an annihilation operator. We define 
the matrix elements CM/^] as 



V{r)cl = ^2c i /3 C[l3\a]ri(f}, 



(12) 



where a and (3 specify a one-body electronic state. For 
sake of simplicity, we will make the r dependence of CM/3] 
implicit from now on. Since rj(r) is nonlocal, it does not 
commute even with creation and annihilation operators 
at different positions (r'). From the expression of 77(f) in 
eq. ([9]), we can readily see this, 



rj{f)c* (r') = ((x' — x, y' — y)c^ (r')rj(r). 



(13) 



On comparing this with eq. (fT2")) . we also see that C[a|, 
is diagonal in this real-space representation, 



C[r' 1 \r' 2 ] = S(r[-r' 2 )ax' 1 -x,y' 1 -y). 



(14) 



In the momentum-space representation, however, a and 
(3 is specified by the crystal momentum k and the band 
index n — 1, 2, • • • , Nb- The creation operator of such a 
Bloch state is defined as follows, 



n.k 



1 

N 



X,fc(^)4> 



where u n g(r) is the periodic part of the Bloch wavefunc- 
tion. Then, substituting the above equation into eq. (|13p , 
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we readily obtain the explicit matrix elements for CM 
in this momentum representation, 



C[n, k\n\ k'] = (i- l)8 n<n 'S^ g 



+e -mv+vo)^ Un ^ kv) 



(15) 



where the inner product between \u n g) represents an in- 
tegral over the unit cell. These inner products are noth- 
ing but the gauge connections in k space. Note that, 
in this momentum representation, CM/3] is no longer di- 
agonal. Instead, it takes a matrix form representing a 
k-derivative or k-covariant derivative. The above rela- 
tions defined in eqs. (fT2|) -(fT5 |) play a fundamental role in 
the reminder of the present paper. 



1. The determinant formulae 

Using the ground state wavefunction defined in 
eq. (fTTj) . let us formulate a general expression for its ex- 
pectation value of the DOP, 



(r?(f)) = (GS|»?(f)|GS>. 



(16) 



Noticing that for an arbitrary permutation of N e quan- 
tum numbers, 



sgn(P) = 

(0|ciC 2 • • ■ CN e -l c Ne C p(N e ) C P(N e -l) 



L P(2) L P(1) 



i)|0), 



one can readily verify, 

( V (r)) = (0|77(r)|0)^sgn(F) 
p 

xC[P(JV )|iV e ]-..-C[P(2)|2].C[P(l)|l], (17) 

where (0|?7(r)|0) is a contribution from the uniform back- 
ground p. Leaving further analysis on (0|?7(r)|0) to the 
following subsection, let us for the moment consider 
eq. (fT?]) . The summation over P should be taken for all 
the possible permutations of N e occupied states. Here, 
let us use eq. (fl5j) as CM/^]> bearing in mind a band 
metal/insulator. Then, in order to interpret the summa- 
tion over P as a determinant, we will arrange or rearrange 
every row and column in such a way that all the matrix 
elements CM/3] with a and f3 being occupied Bloch states 
should be accommodated in the upper-left block of C- Let 
us call this N e x N e submatrix as C', i-e., 



c = 


f 

c 


\ 

* 









Then one can reinterpret eq. (fl7| as 

(V(r)) = (0|r?(r)|0) detC'. 



(18) 



Note that the final formula (fl8|) itself is independent of 
the representation. In fact, when |GS) were given by a 
Slater determinant composed by atomic orbitals totally 
localized in the real space, then we would use eq. (fl4|) 
instead of eq. (fTS"|) as C and thus C'- However, a non- 
interacting band metal/insulator wavefunction (w.f.), in- 
cluding gapped mean- field w.f., is usually composed by 
extended Bloch w.f.. Thus, we will evaluate eq. (|18|). 
using mainly its momentum representation. 



2. Contribution from the uniform background 

Before ending this section, let us estimate the contribu- 
tion from the uniform background, i.e. (0|?7(r)|0). From 
the definition of eq. ([9|), it reads, 



log(0|77(r)|0) 



-p^2 lo gc( x ' 

x' , y' 



x,y'-y). (19) 



By using the real-space representation of the matrix C, 
i.e. eq. (JHJ), it could be rewritten as 

log(0|77(r)|0) = -plogdetC 



Recall that p is the average particle number per site; 

P = N b N' J 

v 



is. in terms of the filling fraction per site 



|?7(r)|0) has the following compact form, 
1 



(0\ V (r)\0) = 



(detC)"' 



(20) 



We will use this expression in section. IV. 

Apart from its compact re-expression, eq. (|19|) itself 
can be directly evaluated. Notice that we have only to 
estimate the amplitude of (r)(r)}, so as to check whether 
it (does not) vanishes in a metal (insulator) or not. Thus 
we will focus only on the real part of log(0|?7(r)|0). At 
the leading order in the thermodynamic limit, i.e., at the 
order of 0(L 2 ), the summation with respect to (w.r.t.) 
r' = (x',1/ 1 ) can be replaced by an integral (see ap- 
pendix. A), 

iog<o|»,(f)|o> = -^ 



27T /■27T 







(I 



dO x d9 y log [ {e ie * +i) + (e w « - 1) 



jo* 



and z' 



In term of complex variables, i.e. z = e" 
e l ® v , this can be further written into the double contour 
integral around a unit circle, 



log<0|rj(f)|0) 
= -N, 



dz 

1*1=1 27r * 



dz' log [z' - zi(z)] 



\z'\ 



= 1 27TI 



zz 



(21) 
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In the r.h.s., we introduced z\(z) = — z + 1 — i. When 
|zi(z)| > 1, the integrand of eq. (j2"Tj) . seen as a function 
of z' , has no more poles other than z' = 0. Thus z' 
can be trivially integrated away. In other words, if we 
decompose the contour |z| = 1 into C> and C<, where 
C> = {z\ \z\ = l,\zx(z)\ > 1} and C< = {z\ \z\ = 
1, |^i(2)| < 1}, then one can evaluate the integral along 
C> in the following way; 



f — I 



dz' log [z' - Zi(z)] 



l c> 2ni 7j a /| =1 2ni zz' 
dz log [- zi{z)] 



2ni 



z 



= « lo g 2 +-E 



sin(n7r/4) 9 

1 2 L _j_ j 7^ 

7r * — ' n 2 2t 16 

n— 1 



As for the contribution from C<, one can safely verify 
that it is only pure imaginary. We also analyzed the con- 
tributions at the next leading order, i.e. at the order of 
O(L), which also turns out to be pure imaginary (see 
Appendix A). We thus obtain the following estimation 
for the uniform background contributions, 



log I <(%(0|0)| 
r3 log2 



7T 

n=l 

which can be evaluated as, 



1 sin(n7r/4) 



2f n 2 



N e + 0(1), (22) 



\(0\v(r)\0)\ = exp[-O.46JV e + 0(l)]. (23) 

The fact that log |(0|r?(f)|0)| < - l -^N e in the thermo- 
dynamic limit is one of the indispensable ingredients for 
our conjecture to make sense, the meaning of which will 
become clearer in the section. III. 



III. THE BAND METAL CASE 

Our objective here is to demonstrate that the expec- 
tation value of the DOP indeed vanishes for the metallic 
state. By considering a single band metallic case, we will 
see explicitly how the presence of Fermi surface leads to 
the vanishing of the expectation value of the DOP. As a 
consequence, we are led to classify various types of Fermi 
surface (or Fermi sea) into three categories by their topol- 
ogy. In the multiband case, attentions should be paid to 
the role of the gauge connection in the momentum space. 



A. Single band — topology of Fermi sea 



good starting point for understanding why our DOP can 
possibly distinguish between a metal and an insulator. 
The matrix element (fT5j) reduces in the single-band case 
to, 



C[n\n'] = V2e l i n x 

n' ^ n x^' :1 .—l^ny,ri' + Py ^n x ,n' x ^n y ,n' — 1 



(24) 



where 



p(0) _ J_ e iz^ e - i3 c( x + x o) 
= J_ e - i l 7r e - l T L (y+ao) 

y V2 



(25) 



Here we have parameterized crystal momenta as n x r y ) = 
1, 2, • • • , N x ( y y, k x = 2Tm x /L,k y = 2nn v /L. Our Bril- 
louin zone is thus doubly periodic, i.e., n x + N x = 
n x ,n y + N y = n y , forming a torus in fc-space. Although 
N x and N y are same, i.e. N x & — N y a = L, we dare use 
different symbols, N x and N y , only in this section just 
for clarity of the following explanations. The substantial 
simplification seen in eqs. (f2"4")) and (|25p . when compared 
with eq. (fT5|) . is clearly the disappearance of the gauge 
connections. In the single band case, the gauge connec- 
tion associated with an overlap of two Bloch functions 
at different fc-points gives at most a trivial phase factor, 
thus erased by gauge transformations. 



1. Two representations 

In order to give an explicit matrix expression of £ using 
eq. (f2"4"]) , let us introduce two representations. The ques- 
tion is how to order N x N y indices, which we will perform 
in the following two steps: 



(k x , £; y ^representation: The most simple and con- 
venient way to give an explicit expression to eq. (|24|) 
is to order all the fc^'s and fe y 's in the increasing 
order of n x and n y . Since in our matrix represen- 
tation for eq. (|24[) . each row or column corresponds 
to a 2D lattice point (n x ,n y ), we attribute, (i) the 
inner microscopic structure (inside a given block) 
to the index n x or k x , and (ii) the outer block struc- 
ture to the index n y or k y . Then, the explicit ma- 
trix element of £ is given by 



Let us first consider the simplest case, i.e., that of a 
single-band metal with iVj = 1. This case is almost trivial 
in the sense that the gauge connection in the k space 
does not appear. However, as we will see below, it is a 



C = V^e 2 ^ x 

'j(N x N y ) + r (0) E (JV„^) + r (0) s (AT*,JV») 



(26) 
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with being an N x N identity matrix, 



2. Three types of Fermi sea 



and 



/ (Af») 



Si 



(JV X 



s\ 



O W0 \ 



(JVx 



(^Vx) y 



M is an x iV shift matrix defined as 

/ 1 o o\ 

: 1 '■• : 



SI 



(N) 





\1 



(27) 



ically as 



and S 



Y 



•■ 10 

1 
0/ 

can be thus written symbol- 



^(JVx,Af B ) 



The evaluation of the determinant (jT8|) would be trivial 

if all the (iV, l)-components s[ N ^ were filtered in Then 
the matrix £' becomes an upper-triangular matrix and 
its determinant becomes just the product of the diagonal 
elements. The question is whether this (N, l)-element 
survives in £' in the /-representation. This is actually 
dependent on the topology of Fermi sea, which we would 
like to discuss now. 

Let us first consider a Fermi sea of trivial topology, 
which has null winding number along both k x and /co- 
axes (Fig. Ufa)). Recall that a (N x , l)-element of S[ Nx) 
in S^* 3 is the matrix element between (N x ,n y ) and 
(1, n y ). Both of these fc-points are, however, filtered away 
in the /-representation for an arbitrary n y , since our 
Fermi sea does not wind the torus along the fc^-axis at 



all. Thus any (N x , l)-elernents of S[ N "' in 
not survive in In a same way, one can easily see that 
(N y , l)-block of S^' Wh) is totally filtered away and does 
not enter into £' either. As a result, £' becomes an upper- 
triangular matrix, whose determinant is identically 1 up 
to a trivial prefactor, 



do 



detC' = 2-sV 



(28) 



In order to see whether the DOP vanishes or not, we 
have to compare this value with | (0|ry(r)|0) |, i.e., eqs. ((22|) 
and (|23p . We readily obtain, 



log | Mr)) | 



y-v sin(mr/4) 



n 2 2f 



N 



(29) 



Thus we have verified that the DOP vanishes in the ther- 
modynamic limit for the metallic state of trivial Fermi sea 
topology. 



Different subscripts, x and Y, are used so as to 
recall that these two matrices contain a nontrivial 
structure at micro and MA CKO-scopic levels. 

2. Filtered- (fca;, fc y ^representation: In eq. (fl~8|) . we 
gave an expression of the expectation value of the 
DOP in terms of the determinant of N e x N e ma- 
trix This matrix, , is a sub-matrix of Q and 
composed of the matrix elements between occupied 
states. In order to construct we filter each state 
by whether it is occupied or not, i.e., we rearrange 
the order of the states in such a way that occu- 
pied states are in the first N e rows and columns 
in the upper-left block of by keeping the ordering 
w.r.t. n x and n y . We call this representation as /- 
representation in the following. It is apparent that 
these two representations are related to each other 
by an orthogonal transformation. 



Let us now consider a less trivial case, i.e., the case of 
such a Fermi sea as depicted in Figs. Hfb) and (b'). The 
Fermi sea winds the torus either along k x or fc^-axes. 
Suppose that the Fermi sea has a filled strip which round 
the torus along the k x direction as in Fig. Q] (b). Then 
the strip is specified by k y = ^j-n y , where n y takes such 
values as n y = no, no + 1, ■ ■ ■ ,uq + A y — 1 on the strip. 
The width of the strip is therefore determined by A y to 
be 2nA y /L (see in Fig. g[b)). 

The Fermi surface topology along the fcy-axis is, how- 
ever, still trivial. Thus (N v , l)-block of E^'^" 3 is to- 
tally filtered. Then the matrix £' turns out to be block 
upper-triangle in the /-representation. The determinant 
of therefore, can be factorized into the product of the 
determinants of diagonal blocks specified by n y . Fur- 
thermore, for any k y or n y which is out of the strip, the 
"topology" along the A^-axis is trivial. Namely, if one 
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(a) 





(b' ) 



FIG. 4: (a) type-0 Fermi surface with uneven boundary and 
annulus. (b) Type-A or B Fermi surface. The strip, which 
is enclosed by the two dashed straight line, winds the torus 
along fcz-direction. (b') Another type-A or B Fermi surface. 
The presence of annulus inside the Fermi sea does not change 
the type of the Fermi surface. Namely, we have only to replace 
A H in eq. J3T) by A* + A*, (c) Type-AB F.S. with a loop 
which winds around the torus along k x 



considers a row specified by n y in the [k x , k y ) -plane such 
that n y < no or n y > no + A — 1, this row does not 
wind the torus along the k x -axis. Thus all the submatri- 
ces (diagonal blocks) in which are specified by either 
n y < no or n y > no + A y — 1 are always the upper- 
triangular matrix, whose determinant is identically 1 up 
to a trivial factor. Then the calculation of det £' reduces 
to those of A y diagonal submatrices corresponding to the 



r2»fif->). 

ri^S^^] can be readily calculated as 



strip. They all have the form V^e^ (/W 
Since det [7^ 

1 + (-l)^" 1 ^ 03 )^, one finally finds 



detC' = 2— e l 
Observing 



i + (-i)^--i (r(°)) 



N x 



l + { _ l) N x -l (Tx) ^ <1 + 2 .^ <l+ l 

for N x ^> 1, we can readily overestimate |det£'| as fol- 
lows, 



| det C | < 2~ e~. 



(31) 



In the end, we takes the thermodynamic limit, keeping 
the ratio A y /N x fixed (to be smaller than 1). Then 

e«x gives a correction of at most 0(1) in the expo- 
nential of the expectation value of the DOP, | (77(f)) | — 
exp [ — O-llA^g + 0(1)] • Thereby the correction is obvi- 
ously irrelevant and the DOP indeed vanishes again in 
the thermodynamic limit. 



Finally it is also possible that the Fermi sea envelopes 
the torus along both k x and fc^-axes as depicted in 
Fig. Q|c) , where the Fermi sea has four pockets at each 
corner of the Brillouin zone (therefore still metallic). In 
this case the submatrix (f is no longer block upper- 
triangle. It is not impossible to write down formally 
det£', but its estimation needs some numerical analysis. 
We have verified by a simple numerical analysis with sub- 
stantial system size that expectation value of the DOP 
for a band metal having this class of Fermi surface ex- 
hibits an exponential decay w.r.t. the system size. 



B. Multiple band case without filled bands 

When there are more than a single band, the situation 
becomes suddenly complicated, because of (i) the pres- 
ence of gauge connection, and (ii) the inter-band matrix 
elements associated with gauge connection. However, we 
can still give a quantitative argument by simply adding 
Nb— 1 empty bands to the single-band case studied above. 
Namely, in such a situation, we have only to calculate the 
determinant of the matrix elements in the lowest band 
n = 1. In other words, only intra-band (n = 1) gauge 
connection appears in our analysis. Because of this sim- 
plification, we can still use exceptionally the single-band 
equation (j!2o]) , but in eci. (p?o]) should be replaced by 
the following N x N y by N x N y diagonal matrices, 



r (o) 



^)[n|n'] = If) • (« a |« a+ , fe, (32) 



Here the unit vectors in the k space were introduced for 
convenience; e x = ^(1,0) and e y = ^(0, 1). 

Let us first consider a Fermi sea of trivial topology as 
depicted in Fig. [Ha). As we have seen in the previous 
subsection, the submatrix spanned by the occupied 
states, becomes trivially an upper-triangular matrix in 



(30) this case. Thus, there is clearly no room for T 



to 

play a role. The determinant again reduce to 1 up to a 
trivial prefactor 2^ e l ^~ N " . 

The case of second topology such as Figs.0Jb) is prob- 
ably more interesting. The arguments are totally parallel 
until det £' is factorized into A y diagonal block matrices, 
which are now not completely same. They still have a 
similar form as V2e l ^(I (N ^ + T x 0) S[ Nx) ). But, due to 



the replacement of eq. (|32|) , these A y diagonal block ma- 
trices are also modified as, 

where n y takes such values as n y = no, no + 1, • • • , no + 
A y — 1 on the strip. The N x x N x diagonal matrices 

F x N:r \n y ) in the r.h.s. are defined as T x N:c \n y )[n x \n x '] = 

T x Nx ' Nv ^ [n x , n y \n x , n y ]. As the number of the permuta- 
tion associated with the determinants of both matrices 
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in eq. (|33| is same, the determinant of the r.h.s. is also 
readily calculated; 

det[I^+Ti N ^(n y )-S[ N ^} 

N Nx 

= i + (-D"-- 1 (4 o) ) * n <«i,«Ki^).(34) 

n x =\ 

Hence, we have 

detC' = 2^V* JVe x 

n a +A y -l AT X 

n {i+c-D*- 1 ^) 'nwv +? j} 



"x = l 



where n n *=l * s a P r °duct along a closed loop parallel 
to the fc^-axis. One might dub this loop as C(k y ) in 
the sense that it is specified only by k y . The product of 
gauge connection along this closed loop can be rewrit- 
ten in terms of the integral of the abelian gauge field, 



iVx 



no 

n~=l 



exp 



i.Jfel 



l,fc+e T ) 



C(k v ) 



h,k\ d K 



^1 k) 



o 



. i 



(35) 



Since this gauge field is pure imaginary, the amplitude 
of the r.h.s. is always unity in the thermodynamic limit. 
We thus conclude that the estimation in the previous 
subsection, i.e., eq. pip , is still unchanged even in the 
presence of the abelian gauge connections. The DOP still 
vanishes as ~ exp[-0.1L/V e + 0(1)]. 

In the case of the multiple band with finite number 
of completely filled bands, we have no direct analytic 
proof that the expectation value of our DOP vanishes in 
metal. This is because all the filled bands form a torus 
and the situation associated with these filled bands is 
essentially same as the F.S. depicted in Fig. [He). In 
those situations, we need to check our conjecture with a 
help of some numerics in future. Summary of this section 
is listed in the table. II of the section. VI. 



IV. BAND INSULATOR/MEAN-FIELD 
ORDERED STATE 

Observing the results in Sec. Ill that the DOP van- 
ishes in the metallic case whenever it could be evaluated 
so far, we will complete our another task in this sec- 
tion. Our purpose here is to verify that the expectation 
value of the DOP indeed remains finite in band insula- 
tors and gapped mean-field ordered states. In the crystal 
momentum space, the lowest N v bands are completely 
filled. The remaining Nb — N v bands are totally empty. 
Then, on filtering away these empty bands, the expecta- 
tion value of the DOP reads, 



(ry) = (GS|tj|GS) 



det C 
(detC)"' 



Here we used eq. ((20)) as the background contribu- 
tions (0|£|0). C' is the upper left part of £ in the /- 
representation, i.e. spanned only by the filled N v bands. 
Unlike the F.S. depicted in Figs. (Ha), (b) and (b'), it 
was totally impossible, even in N v = 1 case, to factorize 
this determinant into the product of diagonal elements or 
diagonal blocks. Additional complications appear, when 
one considers a more general case, i.e. 1 < N v < Nb- The 
gauge connections between different filled bands make it 
much larger the number of the permutation associated 
with the determinant of 

In order to resolve these difficulties, we first focus on 
the most strongly localized limit, i.e. the atomic limit. 
We define this limit, such that a Hamiltonian is decou- 
pled into a local Hamiltonian defined for each unit cell. 
These local Hamiltonians commute with one another. 
Thus eigenstates of a full Hamiltonian are given by eigen- 
states of each local Hamiltonian, which we call as local- 
ized atomic orbitals. Bloch bands, which is totally flat in 
/c-space in this limit, can be then obtained, on Fourier- 
transforming these localized atomic orbitals. The peri- 
odic part of the Bloch function \u r) thus obtained has 

no k-dependence in this limit. 

Starting from this atomic limit, we will expand (17) 
in powers of the number of k-derivatives of this periodic 
part. This treatment is verified, provided that our system 
is close to the atomic limit. We will perform this expan- 
sion up to the second order and observe that | (77) | in the 
insulating state is characterized by the quantum metric 
tensor defined in eq. (f3"9"| (Sec also eq. (HU)). When inte- 
grated over the crystal momentum, this gauge invariant 
2-form is directly linked to the physical quantity called 
as the localization length £ (eqs. (|4QI4ip ). Then, we dis- 
cuss that our DOP remains finite, as far as this localiza- 
tion length is finite. Before developing this perturbation 
expansion specifically, we give a brief reminder to the 
relation between localization length and metric tensor, 
including their physical meanings. 



A. brief review — arbitrariness of Wannier function, 
its spread functional and quantum metric tensor 

The relation between the quantum metric tensor in the 
k space and localization length were discovered in the 
process of constructing a gauge invariant localized Wan- 
nier function. Taking Fourier transformation of Bloch 
w.f., one can usually introduce a spatially "localized" sin- 
gle particle w.f. called as Wannier functions, 



1 



e ik-[f-R) u 



(36) 



fceBz 



However, this single particle w.f. is clearly not invari- 
ant under the U(l) gauge transformation, u n ^(r) = 

e'*»'« n while any physical quantities should be. 

In the case of a composite set of bands, say N v filled 
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bands, the corresponding gauge transformation is gen- 
eralized into SU(N V ) gauge transformations. In either 
case, shapes of Wannier functions depend on a particular 
gauge fixing and are unphysical by themselves. 

Marzari and Vanderbilt introduced a unique set of 
gauge independent Wannier functions, by minimizing 
the spread functional associated with this single parti- 
cle W.fjii, 



= ]T [(0,n|r 2 |0,n) - ((0, n^O, n) 



n=l 



(37) 



where measuring f2 at R = does not lose any general- 
ity. Namely, by optimizing the SU(N V ) gauge degrees of 
freedom such that this spread functional Q is minimized, 
we could defined the maximally localized set of Wannier 
functions {\R, n)}. This set of Wannier functions is b 
definition gauge invariant. They further prove in ref. 
that, when optimized, this spread functional becomes 



ni = 12 [(o> n i^io>™> - 12 K-K> m M$ 



(38) 



R.m 



The localization length £ is introduced as this minimized 
spread; £ 2 = fti. Using the identity, 

(R,n\r^r v ■ ■ ■ \0,m) = 



-y 



Ak-R 



A 2 



K.fcK jV ^)(* V ^)---| U m.fc)' 



fceBz 



one can then express this length £ in terms of the follow- 
ing "metric tensor" defined a o 12 ' 17 i 19 , 



N v 



9^ 



N v 



12 { 9 ^ U n,i\ U n',%){ U n',%\ dk " U nS) ' ( 39 ) 



n' = l 

Namely, we have 



£2 = V- 1 2 £ 2 = J_ V- 



(40) 



fceBz 



As is clear from this definition, £ vanishes in the atomic 
limit, where \u t) has no fc-dependence. 

The tensor quantity g^ lL , is dubbed as the quantum 
metric tensor, in the sense that it measures the mini- 
mal distance between two sets of N v eigenbasis {\u n g)} 
and {\u n ft +d j~)}- I n fact, when optimizing the SU(N V ) 

gauge degree of freedom such that J2 n =i I ( u n.k\ u n,k+dk) I 
is minimized, this quantum distance is indeed specified 
by g^ v (k)dk^dk^. 

The localization length thus introduced can be general- 
ized into the following correlation function of macroscopic 



electronic polarizations, measured w.r.t. the ground state 
many-body w.f. jGS)^ 2 -, 



£ 2 = ^(GS|(A M -A M ) 2 |GS), 



(41) 



where = Ylr=i denotes the total electron's posi- 
tion operator. This correlation function is infinite when 
a system has a Drude weight or peak, while remains fi- 
nite without any Drude weight^ 2 -. Thus, £ is also finite 
(infinite) in insulators (metal). 



B. Disorder operator from the viewpoint of 
localization length 

We see below that the expectation value of the DOP 
in insulating states is simply related to the localization 
length, as given in eq. ([50)) . To be specific, starting with 
the atomic limit, we expand \og{n) in powers of 9Au t). 
The zero-th order, which corresponds to the atomic limit, 
gives a constant contribution. At the first order, one 
finds a correction due to the gauge field. However, this 
is pure imaginary, and irrelevant to the amplitude of (rj) . 
The first relevant contribution to \(rj)\ appears at the 
second order, which will be interpreted, with the help of 
eqs. (|39l40p . as the square of localization length. Then, 
eq. (|50|) in combination with eq. (j4Tj) tells us that the 
DOP remains finite in insulators, as far as this pertur- 
bation expansion is valid. Since the localization length 
diverges in metals, a naive extrapolation of our second 
order analytic expression is also consistent with our pre- 
vious observations found in the metallic case. 



1. The zeroth order - atomic limit 

In the atomic limit (AL), where \u g) has no fc- 
dependence, the gauge connection becomes a trivial ma- 
trix, i.e. (u z\u , p) — <5 n ,n'- Then the matrix £ takes 
the following decoupled structure in the momentum rep- 
resentation, 

CalK k\n', k'] = S n ,n' ® Csb[*|A']- 

From eq. (|15p. Csb is simply given as follows, 

Csb[£|#] - (7*<%+^,fc/ + + {lyk+e y S> ~ S k,k>)> 

where 7 M is v^e^T^ ' . This A 2 by A 2 matrix is diag- 
onalized by the Fourier transformation, where the unit 
cell index: R — (n x &,n y &), specifies its eigenvector, i.e. 
e %k-R rp ne corresponding eigenvalue of Csb reads, 

\{R) = ( 7 xe ieVi? + i) + (l v e l ^ A - l). (42) 

As the band index and momentum index are decoupled 
in Cal, we can easily filter unfilled bands so as to obtain 
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log det Cal = Tr 1°8 Cal- ^ n other words, we have only to 
sum logA(i?) w.r.t. the unit cell index R and multiply 
the number of filled atomic orbital within each unit cell; 



N 

logdetC AL =^ E 1 °SHR)- 



Notice that X(R) in eq. (|42| is basically same as the orig- 
inal ((R — r) defined in eq. ([7]). Thus we can relate the 
r.h.s. with det Q used in eq. (f2"TJ)) . 

log det Cal = ^ • N b E log A(i?) = v ■ log det C, 



where v 



N, 



N b N 2 



This tells us a rather trivial fact that 



our DOP becomes unit in this atomic limit, (t))al 

det Cal _ 113 
(detC)" 



(0|r/|0) = l/(det£)", when exponentiated. On the other 
hand, using the inverse of Cal' 



Cal l Vi,k\n',k l ]=5n,n'j^^e 



l{k - k '^- R \{R)-\ 



the 1st order term w.r.t. the "gauge field" A^ is given as 
follows, 



N v 



Tr^-V] = E EE4 n|n] w 



E 



N 2 ^ {~j x e iS *- R + (-y y e ig v R - l) 



The coefficient of A^ can be further evaluated by the 
following double counter integral, 



2. First and second orders 

We now introduce weak transfer integrals between the 
localized atomic orbitals, taking into account weak k- 
dependences of \u g) perturbatively. Namely, we de- 
velop a perturbative expansion in powers of the number 
of its fc-derivative. To concretize the perturbative evalu- 
ation of det £'/(det let us consider its logarithm, i.e. 
log det C = TrlogC, and decompose C as C = (' AL +A(' . 
AC describes the deviation from the atomic limit, in 
which the information of weak derealizations are en- 
coded through the "gauge field" A^Qc), 



A^ n \k) = (t 



%'k+ej ^ 



2tt 



= — (u r\dh u , A + 

T \ n,k\ n'.k/ 

5(T) a <«n.sl fl 2.V«> + °(^)-( 43 ) 

AC is written in terms of this "gauge field" as 



A('[n,k\ri,k'] 



where n and n' are restricted within filled bands. Observ- 
ing that AC contains at least the first order of d%\u g), 
we first expand Trlog(C AL + AC) w.r.t. AC, 

Trlog(C AL + AC) - logdet[C AL ] + Tr^^A?] 

- ^[CAL-'ACCir'AC] +Om. (44) 



E 



1 

N 2 4^ (<y x e iS *- R + {j y e iS y R - 1) 



[0,2*] X [0,2*] 



+ i + e' w y - 1 47r 2 



d9 x d9 v „, 1 



L 2J> 



where the complex variable e 18 ^ = 7 M e lef, i? was intro- 
duced. Notice that the integrand in the r.h.s. can be 
simply replaced by —ide^log^ 18 * + i + e l6y — 1]}, 



dO x dOy 



[0,2*] x [0,2*] e ie * +i + e ie « -1 4tt 2 



i 

1^2 



d0 A 



where /I denotes the other coordinate than /i e.g. x = y. 
As was clear from Figs. 12131101 log[e l9x + i + e l9y — 1] has 
a branch cut, which runs from (— ir/2, 0) to (0, — tt/2) in 



the 



plane. Thus a surface term associated with the 



^-integral in the r.h.s. gives 2m, when its integral path 
crosses this branch cut. Then, irrespective of a specific 
choice of the principle branch, this double counter inte- 
gral is always estimated to be 1/4 (see Fig.QIJlin the case 
of /j, = x); 



d9 x dOy 



[0,2*] x [0,2*] e i9 * + i + e™y 



- 1 4tt 2 



i 

'4^2 



dOn ■ 2m' = — . 

F 4 



(45) 



As was shown above, its zero-th order term is can- 
celled by the contribution from the uniform background 



Using this coefficient, the 1st order term w.r.t. A^ in 
eq. (|44[) then becomes as follows, 



12 



fi=x,y 



(46) 



r 



where we further expanded the "gauge field" A^ up to 



0(di), according to eq. (|g3l 

Note that the integrand in eq. ([46|) contains the 2nd 
order derivative of \u n ) w.r.t. the crystal momentum, 
i.e. (u n \dl u n }- As we mentioned above, this 2nd or- 
der derivative term is more important than the 1st order 
one, i.e. (u n \dk^u n ), when it comes to the amplitude of 
(r/). Namely, the latter one is clearly pure imaginary, 
while the real part of the 2nd order derivative term, i.e. 
(dk u n \dk^u n ) , constitutes a metric tensor in combina- 
tion with Tr[ (Cal AC') 2 ], which we will sec below. 

Tr[ (C AL _1 AC') 2 ] has a relatively complicated form: 

^ [Cal ~ 1 AC' Cal _ 1 AC'] = 

N v N v 

EEEE4 1| "' 1 (M n ' |,i, (fc') 

fl % v n— 1 n' — l 



7m £ 



N 2 7 T e lS " A 

R ' 



1 

N 2 7 r e 1 ^'^' 



E 



i + 7 y e ie »'- R — 1 

i(k'+e„-k)-R' 



■R' _ 1 



(47) 



above, 
1 

N 2 



E 



e .(S fl +e, J )-R 



(lx e 



ie x -R 



l + Jy e 

d , 



" ' R — l) 
1 



[0,27r]x[0,27r] 



r 



d9 x d9y 
4tt 2 ' 



Then, integrated by part w.r.t. the r.h.s. clearly 
vanishes when fj, ^ v, while becomes identical to the 
l.h.s. of eq. (|4"5l) when /i = v. Thus, they are indeed 
diagonal w.r.t. [i and v; 



-y- 



e i{e l _ l +e„)-R 



N 



(7x e 



iS^-R . 



lv e 



ie y -R ^ 



1 V / 1 \ 



Using this, the 2nd order term w.r.t. A^ is then esti- 
mated up to the accuracy of 0(d~), 



Tr[c AL 'ac'C; 



I /-I — 1 
AL 



^] = -zE E 



fi—x,y n,n'—l 



dk x dk y | (u^ | fy, u n > | 2 + O ( i , 9f ) .(48) 



However, this can be still simplified substantially, by 
noticing that only the leading term is relevant to our 
O (di) -estimations, in the expansion of A v (k') around 

k' = k, i.e. A u (k') = A v {k) + (# - k)x • d kx A„(k) + ■■■. 
Namely, higher order terms in this expansion lead to at 
least the third order fc-derivative of \u n ^), when substi- 
tuted into Eq. (|47|) . Thus, we fairly replace A u (k') in Eq. 
(l47|) by A v {k). This replacement enables us to sum over 
k! , which brings about a factor N 2 5^ Accordingly, as 
far as the C(92)-estimation is concerned, eq. (|47p reduces 
to the following simple form; 



[Cal _ 1 AC' Cal ^C] 



N v N v 

EEEE4 nK1 W nl « 

ji.v £ n=l n' = l 



7M7.e l( ^ +e "*- ) - fl 



iV 2 



^ {l x e iS *- R + i + 7 y e"V- R - l)' 



Observing this simplification, one can readily see that 
the coefficients of A u ■ A v become diagonal w.r.t. /j and 
v. Evaluate the summation w.r.t. R in the r.h.s., in term 
of the double counter integral w.r.t. 9 X and 9 y introduced 



E 



dkx dky X 



Notice the coincidence between the factor j of eq. (|48|) 
and that of eq. (|4"6"]) . This coincidence actually helps us 
to combine these two contributions so as to obtain the 
gauge invariant ((^-estimation for the amplitude of 
our DOP, 

In | <i7(r)>| = 

N v N v 

E {( d k>u n \d k ^u n ) - \{u n '\d k ^u n )\ 2 y (49) 

n— 1 n' — l 

Namely, the integrand is nothing but diagonal compo- 
nents of the quantum metric tensor defined in eq. (|39|) . 
Using the localization length defined in eq. (f40|) , we then 
rewrite this as follows, 



\( V (f))\=exp[- -~e\. 



(50) 



where £ is the dimensionless localization length, i.e. 
£ = |. As we have already mentioned, the localiza- 
tion length is infinite (finite) when the system has a (no) 
Drude weight or peaki^. Thus, eq. (f5U]) indicates that 
our DOP indeed remains finite in band insulators/gapped 
mean-field ordered states. 
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This 2nd order estimation was, however, derived per- 
turbatively w.r.t. the fc-derivative of the periodic part of 
the Bloch w.f.. Thus eq. (150|) obviously holds true only 
for small £ region closed to the atomic limit. However, 
this simplest relation gives a key to extrapolate our anal- 
yses in this section to the metallic case studied previously. 
Namely, starting from the atomic limit where the local- 
ization length is zero, we slowly turn on the transfer inte- 
gral between atomic orbitals so as to increase £. Eq. ([50]) 
indicates that increase of the localization length reduces 
the amplitude of our DOP. We could further increase the 
inter-atomic transfer such that the band gap eventually 
collapses and a Fermi surface appears. Although this sit- 
uation is clearly beyond the validity of our expansion, the 
indication of eq. (|50|) is still consistent with our analyses 
in Sec. III. Namely, as the localization length diverges on 
the appearance of Fermi surface, eq. ([50]) means that the 
DOP should indeed vanish in such a metallic case. 



V. "TWISTING" BOUNDARY CONDITIONS 
AND THE DISORDER OPERATOR 

To help readers to obtain a transparent interpretation 
of the disorder operator, we will attempt in this sec- 
tion to clarify the similarity between our DOP approach 
and other recipes or criteria for probing insulating be- 
haviors. A standard recipe is to investigate system's re- 
sponse against "twisting" the boundary condition a 10 ' 20 . 
The twisted boundary condition along the /i-direction re- 
quires the many-body wavefunction $(fi, • • • , fj\r e ) obeys 

Jf-ji+I^-^e^f-,?;,-) (51) 

for arbitrary i, where real- valued <f> can be non-zero mod- 
ulo 2tt. With the ground state wavefunction obeying this 
twisted boundary condition, one can take the 2nd order 
derivative of its eigen-energy with respect to 4>. This 
quantity estimated at a finite system size becomes the 
Drude weight, when considered in the thermodynamic 
limit, i.e. L, N e — > oo. 

Notice that this twisted boundary condition can 
be transcribed into the magnetic flux <fi inserted into 
those Hamiltonians obeying the periodic boundary condi- 
tion 10 ^ 20 . Namely, as far as its eigen-energy is concerned, 
we may as well consider the Hamiltonian defined on the 
torus having its genus pierced by the non-zero flux </>. 
In the following, highlighting the Aharonov-Bohm phase 
created by the DOP, we will argue that applying the DOP 
"successively" along a certain path E on the 2D plane is 
also equivalent to inserting 2ti x (integer) flux into a sys- 
tem; 

JJ r,{r)~e l ^ m ^ r >^\ (52) 

re path E 

Here the r.h.s. is the unitary transformation, sometimes 
dubbed as "twist operator" , which introduces 27r x m flux 
through the genus associated with the /z-direction. 



A. flux insertion and the disorder operator 

To determine the path E in eq. (|52[) and also /i in its 
r.h.s., let us recall the detailed analytic properties of the 
DOP first, 

77(f) = exp logC(f -f){/?(r') -p}] 

x',y' 

({? -r) = -i^TO'-z) - 1) + (c't (v'-y) - 1). 

Here we have already chosen as io = Vo = 0, for the 
sake of simplicity. Its essential feature is actually the 
existence of vortex-antivortex pair. The vortex is lo- 
cated at (x',y') = (x,y), whereas the antivortex at 
(x', y') = (x + L/4, y — L/4). A branch cut runs, e.g., on 
the straight line, x' + y' — x + y, and connects the vortex 
with the antivortex, encoding the U(l) phase holonomy 
associated with those singularities. In fact, one can see 
from Eq. (fT3| that ImlogC(r' — f) corresponds to the 
"AB phase" which an electron at f subjected under 77(f) 
acquires while traveling on the 2D plane. For example, 
under the influence of the DOP, an electron winding the 
vortex anti-clockwise acquires 2n AB phase. Based on 
this observation, we naturally come up with the idea of 
representing the analytic structure of 77(f), i.e. a pairwise 
form of the vortex and anti- vortex, by introducing a flux 
tube piercing the system at these singularities. 

In order to materialize this point, imagine that our 2D 
periodic system, forming a torus, resides in an 3D space. 
Our flux tube also lives in this space and intersects with 
the 2D plane at these two singular points. Namely, on the 
(x 1 , y')-plane, it appears outside the torus at (x,y), i.e., 
at the vortex, and after traveling above our 2D plane, it 
disappears inside the torus at the antivortex (2; + L/4, x— 
L/4), forming eventually a closed loop behind the system 
(see Fig. [Ha) . The viewpoint of Fig. [5] is located clearly 
outside the torus.) 

Let us now consider a successive application of 77(f) 
on the straight line E parallel to x + y — constant, i.e. 
Y\p eE ri(r). In the language of this flux tube, this cor- 
responds to a superposition of the flux tubes which are 
shifted by one lattice site (diagonally along E) between 
one another (Fig. Efb)). When these tubes are super- 
posed, cancellations between vortices and antivortices oc- 
cur in the vicinity of 2D plane, and the singularities on 
the plane totally disappear. As a result, we are left with 
two decoupled vortex lines: one running above, the other 
beneath the system along E (Fig. [|3a)). These are a 
bundle of N/4 quantized (2tt) flux tubes; N/4, because 
of N sites (operators) along E and each operator having a 
branch cut of length, equal to 1/4 of the whole trajectory 
of E (= V2L). 

These two vortex lines can be deformed freely out- 
side or inside the torus without intersecting the torus. 
Namely, as for an electron living in the 2D plane (torus), 
this freedom itself corresponds to a trivial U(l) gauge 
degrees of freedom. Using this freedom, one can readily 
see that these two decoupled vortex lines indeed com- 
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(a) 




Cancellation -471 flux bundle 



(b) 

FIG. 5: (Color online) (a) A quantized (2tt) flux tube asso- 
ciated with the disorder operator (DOP) 77(f) is shown as a 
red(thin) arrows on the (x' , y')-plane. Corresponding to the 
phase inconsistency of log £(r' — f), the flux tube appears out- 
side the torus at (x,y), i.e., at the vortex, and after traveling 
above the branch cut, it disappears inside the torus at the an- 
tivortex, (x + L/4, x — L/4), forming a loop, (b) The path E is 
a straight blue (dotted) line along which a series of flux tubes 
associated with the DOP 77(f) is considered. Flux tubes asso- 
ciated with 77(f) at each lattice site F on E are superposed. In 
this summation over ron E, cancellations between flux tubes 
heading outside and inside the system occur on the 2D plane 
(shadowed region). As a result, those flux which penetrate 
the system completely disappear, leaving two bundles of flux 
tubes above and beneath the system in parallel with E (See 
also Fig. [6ja)). 




271N/4 
(b) 

FIG. 6: (Color online) (a) Two flux bundles and the periodic 
2D flat surface (= our system). The red(solid) line represents 
a bundle of N/4 quantized flux tubes which flows above (<-> 
outside the torus) and in parallel with our 2D plane. The 
red(broken) line is the other bundle which flows beneath (<-> 
inside the torus) and in parallel with the plane, (b) Two flux 
bundles and the torus (= our system). The former pierces the 
handle of the torus encircled by the j/'-axis, whereas the latter 
penetrates the interior of the torus wound by the a/-axis. 

pose a pair of "knotted" rings, consisting of a "Hopf link" 
(Fig. [HJb)). In other words, we have 

W). (53) 

This formula is almost correct, but there are two mi- 
nor imperfections. First, the two sides of eq. (|53p can- 
not be completely identical, since our DOP is not uni- 
tary, whereas the r.h.s. is clearly unitary. We, there- 
fore, slightly soften our statement in such a way that the 
unitary part of fifes 'j(^) ^ s equivalent to the r.h.s. in 
Eq. (|53|) . Secondly, since we freely deformed two vortex 
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lines without intersecting our 2D plane, we can not fix 
in eq. (|55|) a trivial U(l) phase factor e LtpE with a real- 
valued continuous periodic function <p E (r'); 

1. d x dyip E = d y d x tp E , 

2. <p E {x' + L, y') = ip E (x', y' + L) = (p E (x', y'). 

Then, to the best of accuracy, our claim can be reformu- 
lated, in the language of the phase part of the DOP and 
twist operators, as 

J2 ImlogC(^ - r) = y X — {x 1 + y') + <p B {?). (54) 








A 





A x' — y' plane, where branch cut of 
r) runs from (x,y) to {x + k,y—k). The shaded 
while the other side is for — (Right) : 



FIG. 7: (Left) 
Imlog £(r 
side of it is for 

One can easily convince oneself that eq. (|57|l is independent 
from a specific choice of the branch cut. 



B. Direct calculation of phase holonomy 

Eq. (f5"4"]) elucidated by the heuristic argument can be 
directly verified, on taking lattice constant "a" infinitesi- 
mally small with Na, = L fixed. In this continuum limit, 
the summation over r — (x,y) along E is replaced by the 
integral w.r.t. I = with fixed 1= 

dl Imlog^r' — r) 

(55) 



E 



N 



Imlog ^{r 1 — r) 



.2 



V2L„ ^ 
F(x',y';l), 



which is a function of I' — I in general. Then, by in- 
tegrating d r '^F along the redirection over [— -j, -j], we 
can directly calculate the AB phase which an electron 
acquires each time it travels around the periodic system 
along r'^-direction. Let us see this phase holonomy, tak- 



ing to be 



first, 



OF 



dx' = 



■J2 



dl 



d_ 

dx 



yjlmlogC^-r)} da/. (56) 



The integral in the r.h.s. clearly vanishes due to the 
periodicity of log£, in the case of y' — y [— -§,0] (see 
the path B in Fig. [7]). However, as for y'-ye [— -j,0], 
the path of the integral w.r.t. x' (the path A in Fig. [7]) 
inevitably crosses the branch cut where km log £ jumps 
by 2tt: Imlog £(r" — r) has a branch cut emitting from 
f = (x,y) to r* = (x + j,y — j). Thus, in the latter 
case, the surface term gives 2iri, on integrating w.r.t. x' 
over [-f i.e. 



As a result, eq. (|56j) reduces to 



^{nr-D}*-- 



V2L 



-(V2y'-l) 



dl 



,9?r 



(v^y'-O-^ 



2ir. (57) 



This unambiguously dictates that F(V — I) differs from 

8 <p"(x' , y'; I) whi 
<p"(x' + L,y';l) 



2%r- only by those continuous functions ip"(x' ', y'; I) which 



are periodic along the x'-axis, i.e 
<p"(x>,y';l): 



F(l>-l) = ^+<p"(x',y>-J). 



(58) 



On repeating the same argument for y', we can easily 
find that these periodic functions ip" is in turns different 
from -22- only by those continuous functions (p'(x',y'; I) 
which are double-periodic both along the x-axis and y- 
axis: tp'{x' + L,y';l) = ip' (x' , y' + L;l) = (p'{x' ,y';l). 
Then, we safely verified the continuum limit of eq. (154 



dl Imlog — r) 



2tt(x' + y') 
L 



+ if'(x' + y';l). 



(59) 



Let us remark that the path dependence of the r.h.s., 
i.e. I = ^^-^-dependence, only enters into the double- 
periodic function ip' and does not change the strength of 
the flux inserted along x and y-direction. To summarize 
this section, we argue that the successive application of 
the disorder operator (DOP) along those arbitrary paths 
which parallel to its branch cut is equivalent to twisting 
the boundary condition of the system. 



when y' 
when y' 



ye [- 
vt[- 



f,o] 

f,0] 



VI. SUMMARY AND DISCUSSIONS 

Based on the duality picture between 2D supercon- 
ductivity and insulator, we proposed in this paper the 
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disorder operator as a possible candidate for an order 
parameter of the 2D insulating state, which is valid for 
arbitrary 2D lattice models. Namely, the change of phase 
of our disorder operator along a closed loop is nothing but 
the doped particle density inside this loop. Thus, one can 
naturally introduce the disorder operator (DOP) as the 
dual counterpart of the superconducting order parame- 
ter. 

To inspect the validity of this conjecture, we evaluated 
the expectation value of this non-local operator in band 
metal (Sec. Ill) and in band insulators/gapped mean- 
field ordered state (Sec. IV). Thereby, we observed that 
the expectation value of the DOP actually vanishes for a 
wide class of band metals in the thermodynamic limit. As 
for band insulators, we estimated its expectation value 
perturbatively w.r.t. a weak fc-dependence of periodic 
parts of Bloch w.f., bearing in mind the atomic limit 
or weak delocalized insulating state close to this limit. 
Thereby, we found that the DOP is characterized by the 

,2-2 

localization length £ as e~~^ . This localization length 
is always finite in insulators, while diverges in the pres- 
ence of F.S. Thus, these two theoretical observations, i.e. 
in the metallic case and in the weakly delocalized insu- 
lating case, naturally lead us to speculate that the ex- 
pectation value of the DOP in insulators could be always 
characterized by the same form, although our derivation 
in the insulator side is valid only for small £. This specu- 
lation must be checked in future, with a help of numerics 



or another analytic scheme of evaluating determinants. 

One might also wonder about the behaviour of the 
DOP in other electronic states, such as superconductors 
(SC). As for these off-diagonal long range order states 
(ODLRO), we expect that its expectation value should 
also vanish^£&£&^, from the simple quantum mechan- 
ical commutation relation between a local density and 
creation operator, i.e. [pf^ctj = ^/jct . Namely, the 
particle number at each site pp and the phase part of the 
SC order parameter, e.g. c^ct ^, are canonical conjugate 

to each other. Thus, in such an ODLRO where ct^ct^ is 
condensed uniformly, the particle number for each site is 
totally undermined in an entire system. Since the local 
electronic polarization a(r) and phase part of the DOP, 
i.e. lnry(r), are given in terms of the linear combination 
of these local electron densities, they should be also in- 
definite. In other word, the DOP should not acquire a 
finite amplitude in these ODLRO states. 

The standard approach of detecting insulating states 
(and also superconducting states) is to measuring the 
ground state energy variation w.r.t. the infinitesimaHy 
small Aharonov-Bohm (AB) phase inserted in parallel 
with system o 10 ' 18 ' 19 ' 20 . In order to make a connection to 
these conventional approaches, we prove that applying 
the DOP successively along a closed path also eventually 
ends up with the AB phase inserted in parallel to the 2D 
plane. 



ground state w.f. 


band metal 


band insulator/gapped m-f state 


type of F.S./ localization length 


type. O 


type. A or type. B 


type. AB 


general £ 


small £ 


single band 


e -o.u.JV e 










m-b. w/o. filled bands 


e -0.11-iV e 


< e a • (.-a-"-*. 






-zip 


m-b. w. filled bands 













TABLE II: Summary of Sec. Ill and IV. "m-b." denotes "multiple band". Type. O, type. A(B) and type. AB denotes a type of 
the F.S. of 2D metal (See Fig.©. £ = * denot es a localization length along the /^-direction (see eqs. (|40I41|I for its definition). 
"— " is the case which we cannot access by our analytical scheme. 



Observing these several circumstantial evidences in- 
cluding the connection with the conventional approach, 
we believe, in spite of its complexness stemming from its 
non-locality, that the DOP is finite only in an insula- 
tor. Then, turning back to our original motivations, i.e. 
microscopic identification of the counterpart of magnetic 
penetration depth and coherence length, we are now al- 
lowed to push forward a naive thought on this primary 
motivations. We expect that the localization length £ 
characterizes not only the expectation value of the DOP 
in insulating states, but it also specifies the counterpart 
of the magnetic penetration depth, by the following rea- 
sons. Notice, from eq. (PETj) , that the localization length 



measures how easily the localized electrons constituting 
the background non-doped insulating w.f. could be po- 
larized, when an external electric field is applied, or in 
other words when a test charge is introduced into this 
system. When this test charge is regarded as a sin- 
gle (few) doped particle, the polarizability of the back- 
ground insulating electrons then could be transcribed 
into how easily a single doped particle push away these 
background electrons and acquires its own seat within a 
bulk, or in other words, forming a vortex within a bulk. 
Thus, we expect that there should be a certain amount 
of positive correlation between the localization length £ 
measured w.r.t. non-doped gapped mean-field ordered 
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w.f. and the counterpart of the magnetic penetration 
depth in its weakly doped regions. These speculation 
could be directly judged, when one construct a Ginzburg- 
Landau type theory for doped 2D insulating state, by 
using this non-local field, or alternatively, when calcu- 
lating its multi-point correlation functions, w.r.t. the 
non-doped insulating w.f.. 
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APPENDIX A: CONTRIBUTION FROM THE 
UNIFORM BACKGROUND 

When we prove our DOP indeed vanishes in a band 
metal, it is crucial that the contribution from the uni- 
form background (0|?7(r)|0) overkills the contribution 
from detC', i.e. |(0|r?(r)|0)| < e^^^". Therefore, in 
this appendix, we will estimate the contribution from the 
uniform background, 



f' 

ln[ - l ( e l ^ x '~ x - x ^ - 1) + ( e l ^ y '- y - ya) - 1)], 

(Al) 

up to the order of 0(1) and prove that this is indeed the 
case. This summation can be replaced by the integral, as 
far as 0(L 2 ) is concerned: 



ln(0|?7(f)|0) 



Nn 



N b £ 



N b N 2 

ll ,y — a, 

ln[ - i( e l ^ {x '- Xo) - 1) + (e l ^ {v '- ya} - 1)] 

/ d6 x / dO y In[t + e ie * - 1 + e ie «] 
Jo Jo 

dz f dz' bx[z' — zi(z 
| z |=i 2tti J\ z '\=i 2th z ■ z' 



4tt 2 
N,, 



, (A2) 



where we introduced complex variables z = e % ® x , z' = 
e y and z\{z) = —z + 1 — i. The additional factor Nf, in 
the 1st line comes from the summation within the unit 
cell. Then, we will first divide the above z-integral into 
the following two parts, according as the position of 21(2) 
in the complex plane; 



Namely, the region G> i < are defined as follows (see 
Fig.©: 

C< = {z I \z\ = 1, \z x {z)\ < 1}, 
C> = {z I |z| = l, \ Zl (z)\>l}. 



We can easily see that the 2nd term in eq. (|A3j) does 
not contribute to the real part of eq. (|A2j) , by choosing 
the branch cut of the logarithm, i.e. log[z' — zi], so that 
it depends on z%. Namely, let it run from Z\ to infinity 
passing through t^t (see in Fig. O . Note that the real 

part of eq. (|A2|) is free from the specific choice of the 
branch cut. According as this choice, we will decompose 
the z'-integral into the following four parts: 

where the closed path C is the bold loop depicted in 
Fig. [Hand 

Ci = {z' I arglnz' = arglnzi + 0, \z x \ < \z'\ < l}, 
C 2 = \z' I \z / -z 1 \ = +e}, 

C 3 = {z' I arglnz' = arglnzi - 0, |zi| < \z'\ < l}. 



First of all, the third term in the l.h.s. of eq. (|A4[) clearly 
vanishes, when e is taken to be infinitesimally small: 



dz' \n[z' - Zl ] i e f 27r 



C-2 



2-ki 



2tt zi Jo 



e l U, (A5) 



where z' — Z\ = ee 1 ^. As the integrand in eq. (| A2|) . seen 
as a function of z', has no more pole than z' = inside 
C, the first term in eq. (IA4j) reads: 
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FIG. 9: Dotted line emitting from Z\ to infinity is the branch 
cut of ln[z' — Zi] which we have chosen in eq. (|A4[) . arg ln[z' — 
Zi] takes +0i in the shaded side of this dotted line while 2n 
in the other side. 




FIG. 10: (Left): A 9 x -9 y plane, where branch cut of ln[i + 
e ie * -l + e i£ *»] runs from (-§ , 0) to (0, -§). The shaded side 
of it is for —j— while the other side is fo r — + . (Right): 
One can easily convince oneself that eqs. (|45IA12p are inde- 
pendent from a specific choice of the branch cut. 



Due to our choice of the branch cut, the second and 
fourth terms in eq. (|A4| can be substantially simplified: 



dz' m\z' -Zi\ + 0-i 



2m 



dz' ln\z' - zi\ + 2m 



c 3 



2m 



dz'— = lnj^ij 



(A7) 



Then, by combining eq. (|A6|) and eq. (|A7|1 . we finally 
obtain the 2nd term of eq. (|A3[) as, 



dz 1 



dz' ln[z' — zi] 



c< 2m z Jj 2 /| = i 2m z' 



1 



dz 



dk = ^-i. 

2 J c< z 2 7_ a 4 



(A8) 



Thus, as far as the real part is concerned, there is no 
contribution to eq. (|A2j) . when z is located on C<. Then 
we will concentrate on the case where z is located on C> . 
When z on C>, z'-integral can be easily done; 



dz 1 



dz' ln[z' - Zl (z)] 



c> 27ri z 7j z /|=i 2?ri 



dz ln[- l + i + z] 



2?ri 



3, „ 1 v 1 , 1 s n . ,nir, ,9n 
-ln2 + -> — ( — ) sm — )+i — 
8 7T ^ n 2 K J2> V 4 ; 16 



0.464847699170805... + i^. 
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(A9) 



Then, as far as the order of 0(L 2 ) is concerned, logarithm 
of | (0 1 rj(r) 1 0)| can be estimated to be 

ln|(0|?7(r)|0)| = -0.464847... x N e + 0(L). (A10) 

Next, we will estimate the 0(L) contribution to 
ln(0|^(r)|0), which turns out to be also pure imaginary. 
O(L) contribution comes from the error due to the re- 
placement of the summation given in eq. (|A1|) by the in- 
tegral given in eq. (|A2|) . This error can be then estimated 



by the summation of the first derivative of integrand in 
eq. (|A2[) w.r.t. 8 X and/or 6 y . In other words, its leading 
order correction is always proportional to the following 
integrals with some real-valued coefficients, 



2-7T /■2-n- 



d0 x d0 y d 
4tt 2 d6 u 



+ e l 



9 »]},(A11) 



where fi = x,y. 

Taking the branch cut of ln[i + e l9x — 1 + e l8y ] as 
in Fig. 1 101 one can easily understand these correction 
are pure imaginary. Let us see this, taking /i — x for 
example. In the case of 9 y [— f , 0] (see the line B 
in Fig. ITU|) . the ^-integral always vanishes due to the 
periodicity of ln[z + e 10 * — 1 + e l6y ] . On the other hand, 
in the case of 9 y E [— f,0], the integral w.r.t. 9 X meets 
27r jump when its integral path crosses the branch cut 
which extends from (—^,0) to (0, —5) (see the line A in 
Fig. [TOf . Then, when integrating w.r.t. 9 X over [— 7r, w], 
we obtain 2m as the surface term in the latter case. As 
a result, the above integrals with \i — x reads: 



1 

4^2 



d9 v x 2ni = -. 



(A12) 



In a same way, we can easily see that the integral with 
fj, = y also turns out to be j. After all, the error which 
stems from replacing the summation of eq. (|A1[) by the 
integral given in eq. (|A2|) can have a contribution of the 
order C(L), but it is always pure imaginary. Thus, when 
it comes to ln|(0|?7(r)|0)|, there is no 0(L) contributions. 

To summarize this appendix, ln|(0|?7(r)|0)| is estimated 
up to the order of 0(1) and it reads: 

ln|(0|?7(r)|0)| = -0.464847... x N e + 0(1), (A13) 

which is actually less than — ^^N e . 
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